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We investigate the structure of vortex-bound states in spin-singlet chiral superconductors with (d^ 2 _y 2 zb 
idxyfwsLVQ and {dxz ± idyz)-wsLyQ pairing symmetries. It is found that vortices in the {dxz ± idyz)-wsLyQ state 
bind zero-energy states which are dispersionless along the vortex line, forming a doubly degenerate Majorana 
flat band. Vortex-bound states of {dx2_y2 zb idxyfwsLVQ superconductors, on the other hand, exist only at flnite 
energy. Using exact diagonalization and analytical solutions of tight-binding Bogoliubov-de Gennes Hamiltoni¬ 
ans, we compute the energy spectrum of the vortex-bound states and the local density of states around the vortex 
and antivortex cores. We And that the tunneling conductance peak of the vortex is considerably broader than 
that of the antivortex. This difference can be used as a direct signature of the chiral order parameter symmetry. 


I. INTRODUCTION 


Chiral superconductors are attracting growing interest be¬ 
cause of their potential use for novel superconducting de¬ 
vices and quantum information technology. These uncon¬ 
ventional superconductors exhibit pairing gaps whose phase 
winds around the Fermi surface in multiples of 27r, leading 
to a non-trivial wave function topology and a breaking of 
time-reversal symmetry. The non-trivial topology gives rise 
to a multitude of interesting phenomenair-, in particular sub¬ 
gap states in vortex cores and protected gapless edge modes 
that can carry quantized thermal current and particle current. 
Probably the most prominent example of a chiral supercon¬ 
ductor is the spin-triplet (px ± ip^)-wave state, which is be¬ 
lieved to be realized in Sr2Ru04i24i2, in the A phase of super¬ 
fluid and in two-dimensional cold atomic gasesi^. 

Spin-polarized (px =b ip^)-wave superconductors support non¬ 
degenerate Majorana zero-energy modes localized at vortex 
coresi^“— . These Majorana quasiparticles obey non-Abelian 
statistics and can therefore be employed to implement topo¬ 
logical quantum computing^LS^. 

Another example of a chiral superconductor is the spin- 
singlet chiral d-wave stated. The non-trivial topology of this 
phase is analogous to that of the chiral (px =b ip^)-wave su¬ 
perconductor. However, due to the conservation of spin- 
rotation symmetry, the edge modes of spin-singlet chiral su¬ 
perconductors carry besides a thermal current also a well- 
deflned quantized spin current. Recently, it has been pro¬ 
posed that graphene doped to the van Hove Ailing is a po¬ 
tential experimental realization of the spin-singlet chiral su¬ 
perconductor—. Other candidate materials for spin-singlet su¬ 
perconductivity with broken time-reversal symmetry include 
SrPtAs2i~i2i, the heavy fermion system URu2Si22^“— , and Cu- 
doped TiSe 2 ^~— . 

In this paper, we investigate the energy spectrum and the 
wave function profile of vortex-bound states in spin-singlet 
chiral superconductors with (dx2-y2 zb id^^^j-wave and (dxz i 
idyz)-wsiYC pairing symmetries. Interestingly, we And that 
vortices in the (dxz i idyz)-Y/ 2 iYe state support zero-energy 
states with a flat dispersion along the vortex line (Fig.[T]). The 
(dx2-y2 zb idxy)-Y/SiyQ state, on the other hand, supports vor¬ 
tex bound states only at flnite energy. We show that for both 


pairing symmetries the tunneling conductance peak of the vor¬ 
tex is about twice as broad as that of the antivortex (Figs. [5] 
and[6l). This property is present even at temperatures consider¬ 
ably higher than the energy spacing between the vortex-bound 
states, and can be used as a direct probe of time-reversal sym¬ 
metry breaking and chiral order parameter symmetry. 

The remainder of the paper is structured as follows. In 
Sec.iniwe introduce the Bogoliubov-de Gennes (BdG) Hamil¬ 
tonian of the spin-singlet chiral superconductors in the pres¬ 
ence of a vortex/antivortex pair. The analytical solutions of the 
vortex-bound state wave functions are derived in Sec. lIII Al In 
Sec. lIIIBI we present the numerical results for the local density 
of states around the vortex and antivortex cores and discuss the 
asymmetry between vortex and antivortex bound states. Our 
conclusions and a discussion of implications for experiments 
are given in Sec. [IVl Some technical details of the derivation 
of the vortex-bound states are provided in Appendix lAl 


II. BOGOLIUBOV-DE GENNES THEORY 


At a phenomenological level chiral d-wave superconduc¬ 
tors can be described by the 2 x 2 BdG Hamiltonian V, = 




Hu 


V^k ~^-k/ 


( 1 ) 


and the Nambu spinor <l>k = Here, (Cj^^) 

represents the electron creation (annihilation) operator with 
momentum k and spin s. The normal state h\^ = t(coskx + 
cos ky) tz cos kz — fi describes electrons hopping between 
nearest neighbor sites of a tetragonal lattice, where t and tz 
denote the hopping integrals in the xy plane and along the 
z axis, respectively, and p is the chemical potential. In the 
following we focus on quasi-two-dimensional systems with 
tz ^ t and consider two different spin-singlet chiral paired 
states, namely the (dx 2 -y 2 zb idxy)-Y^di\Q state described by 


Ak = Ao (cos kz + 4) (cos kx — cos ky fzi sin kx sin ky) 

(2a) 

and the (dxz i idy;s)-wave state give by 

Ak = Ao (sin kx sin kz fzi sin ky sin kz ), (2b) 
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where Aq denotes the superconducting gap energy. The su¬ 
perconducting order parameter for both pairing symmetries 
exhibits point nodes at the north and south poles of the Fermi 
spheroid. The gap function (l2bl) has in addition a line node at 
the equator of the Fermi surface, see Fig. [T] The point nodes 
of the (dx 2 _y 2 ± ida:y)-wave state (l2at realize double Weyl 
nodes, whose stability is protected by a Chern number that 
takes on the values ±2.- The low-energy nodal quasiparticles 
near these double Weyl nodes exhibit linear and quadratic dis¬ 
persions along the direction and in the k^ky plane, respec¬ 
tively. This anisotropic dispersion leads to a density of states 
which increases linearly with energy. The point nodes of the 
{dxz i idy;s)-wave state (l2bl) . on the other hand, correspond to 
single Weyl nodes with Chern number ±1—. 

According to the classification of Ref. [sl H\^ belongs to 
symmetry class C, since it satisfies particle-hole symmetry 

CFkC-i = -F_k, (3) 



k /ti 

z 


with C = cFxKL and = —1, but breaks time-reversal sym¬ 
metry. In the following we consider vortex lines along the z 
axis, which allows us to make use of the translation symme¬ 
try along the z direction. Therefore, we can decompose H\^ 
into a family of two-dimensional layers with fixed kz . Hence, 
the analysis of vortex-bound states of the three-dimensional 
superconductor ([T]), reduces to the problem of studying point 
vortices of two-dimensional superconductors as a function of 
kz. By the bulk-defect correspondence of Refs, [sl, S and 
it follows that a two-dimensional Hamiltonian in sym¬ 
metry class C does not exhibit any zero-energy vortex-bound 
states. This is the case for the (dx 2 -y 2 dz i(ia:^)-wave pairing 
state. The superconductor with (dxz i idyz)-wsiYe gap sym¬ 
metry, however, constitutes an exception to this rule. This 
is because for fixed kz Eq. (l2bl) does not have chiral d-wave 
symmetry, rather it exhibits a chiral p-wave character and thus 
belongs to symmetry class D. As a consequence, we find from 
the classifications of Refs. [3, [1, H and [H that vortices in 
the (dxz i idyz)-wsiyQ state support zero-energy bound-states 
protected by a Chern-Simons invariant (see Fig. 


A. Implementation of Vortex/Antivortex Pair 

In the following we discuss how vortex lines along the z 
axis are implemented on a microscopic level. As mentioned 
above, due to translation symmetry along z we can decompose 
the three-dimensional Bogoliubov equations into a family of 
two-dimensional equations parametrized by kz. In order to 
introduce vortex/antivortex pairs we Fourier transform the kx 
and ky components of Eq. o into real space. This yields a 
two-dimensional lattice Hamiltonian on an A x A square lat¬ 
tice, with coordinates ranging from [—A/2, A/2). In order 
to suppress possible surface states of the superconductor we 
impose closed boundary conditions in all directions. The vor¬ 
tex/antivortex pair is implemented by applying a radial profile 
h(r) and phase y) to the gap parameter, 

Ao ^ (4a) 


FIG. 1. (Color online) Energy Spectrum of (a) the (d^ 2 _y 2 + idxy)- 
wave and (b) the (dxz + idyz)-wsL\Q pairing state in the presence of 
a vortex/antivortex pair. The vortex-bound states are highlighted in 
red. 


where n is the vorticity of the vortex/antivortex pair. Note 
that the vortex (antivortex) is defined by a positive (negative) 
winding number arg[A(r)] ds about the vortex (antivortex) 
core, where C is a small circle centered at the core. The func¬ 
tion that parametrizes the phase of the vortex/antivortex pair 
is given byil 


0(x, y) = tan ^ 


/ 2ayC(y) \ 
\x^ ^ y‘^ — o? ) ’ 


(4b) 


where a is half the distance between the vortex and the an¬ 
tivortex. To minimize finite-size effects we choose a = A/4, 
such that the distance between the vortex and the antivortex 
is maximized. With this choice, the vortex and antivortex be¬ 
have like isolated vortices for A large enough. The factor 

l/B 


C{y) = 


l + A- 


M 

N 


in Eq. (l4bl) is used to provide con¬ 
tinuity of the gap phase across the closed boundary, while re¬ 
taining the structure of the original gap phase around the vor¬ 
tex cores. The values A and B are determined by an optimiza¬ 
tion process^. The radial profile of the vortex and antivortex 
at (a, 0) and (—a, 0), respectively, is taken to be 


h{r) = 


0 : 0 < r < 1 

^tanh(r/p) : r > 1 


(4c) 


where p is the size of the (anti)vortex and r the distance from 
the (anti)vortex core. The piecewise nature of the profile is 
used to remove unphysical singularities at the vortex core. We 
observe that the profile h(r) is linear close to the core and 
constant far away from the core. 
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where 



A = ^/A{r){-dl + dl 


2idxdy)^ A{r), 


( 6 ) 


and A(r) = describes an (anti)vortex at the ori¬ 

gin with vorticity n and polar coordinates (r, 0). To simplify 
the analysis of the Bogoliubov equations we rescale the equa¬ 
tions in terms of the characteristic length 

L= 3 , (7) 

which yields the dimensionless variables 


X = ^ and e = emL^ 


e 


( 8 ) 


FIG. 2. (Color online) Comparison between the analytical (solid red) 
and numerical (dashed blue) solutions of the (anti)vortex-bound state 
wavefunctions of the {dx 2 _y 2 -\-idxy)-wsLyQ superconductor, (a), (b) 
Wavefunction amplitude of the lowest energy vortex- and antivortex 
bound state, respectively, (c), (d) Wavefunction amplitude of the sec¬ 
ond lowest energy vortex and antivortex-bound state, respectively. 


III. STRUCTURE OF VORTEX-BOUND STATES 

In this section, we study the structure of the vortex-bound 
states using both analytical and numerical methods. For the 
analytical solutions of the vortex-bound states we focus on the 
{dx 2 _y 2 -f idxy)-v^Siye pairing state. The vortex-bound state 
wavefunctions of the {dxz + idyz)-y/SiYe superconductor can 
be inferred from the published results on vortex-bound states 
of the chiral (px + ipy)-wave stateJ^*^^*^. That is, the vortex- 
bound states of the {dxz + idyz)-wsiYQ stat e (1^ . are obtained 
from the bound-state solutions of Refs. Em an dla^by scal¬ 
ing the gap energy by sm{kz) [i.e., Aq ^ Aq sin(/C; 2 )]. 


A. Analytical solutions for {dx 2 _y 2 -h zda;y)-wave state 


With this, the dimensionless vortex profile becomes 


I A/tanh(I/r/p) 


:0 <r < 1/L 
: r > l/L 


(9) 


For ease of notation, we omit the overbars for the remainder 
of this section. I.e., in the following all variables are assumed 
to be dimensionless. By setting u{y) = exp[i(/ -f ^^^^)0]u{r) 
and i;(r) = exp[i(/ — ^^^)0]v{r), where I = is re¬ 

stricted to half integers (k G N), we obtain for the Bogoliubov 
equations 


--iM+ - 7^ 
-7^m_ + 7^ 


u H —= eu^ 
7 


V H— D_u = ev^ 
7 


(10a) 

(10b) 


with 7 = l/(Aom), M± = |(2 ± 2/ + n), and the second 
order differential operators 


Ls = dl + \dr 


(lOc) 


and 


In order to obtain analytical expressions for the vortex- 
bound states of the {dx 2 -y 2 ida^^)-wave state, we first de¬ 
rive a low-energy continuum description of Hamiltonian (O. 
To this end, we consider a single (anti)vortex at the ori¬ 
gin and assume that the Fermi surface is small, of spherical 
shape, and centered at the T point. Performing a small mo¬ 
mentum expansion, we obtain for the normal state /i(k) = 
^{kl - ky) + - 11 and the gap function^ A(k) = 

/S.<^[{l — k‘l)-\-A]{kl. — ky^2ikxky). Since we consider a quasi- 
two-dimensional system, we can take kz to be a fixed param¬ 
eter and absorb all kz dependent terms in constants. That is, 
we let ji - n and Ao[(l - k^) + 4] Aq. By 

replacing momentum variables by momentum operators, i.e., 
{kx^ ky) —i{dx^ dy), we arrive at the real-space represen¬ 
tation of the continuum Bogoliubov equations 

1,7)0-O' 


D± 


V—p— 


+ 


+ 


h‘^{-l±2l) 


hhA-^ _ 

— 2ft.ft.') dr — h'^dr 


(lOd) 


Analytical solutions to Eqs. (fTOl) can be derived in the limit 
7 1. In this limit the Bogoliubov equations decouple and 

the solutions are given in terms of the Hankel functions of the 
first and second kind, Ha\x) and Ha\x). Thus, we make 
the following ansatz for the wavefunctions 

uir) = fi{r)H^M’^{qr) + 

v{r) = gi{r)H^M\{qr) + (H) 


with q = \/ 27 . The functional form of the coefficients fi{r) 
and gi{r) (with i G {1,2}) is derived in Appendix lAl from 
which it follows that / 2 (^) = /r(^) and ^ 2 (^) = ^i(^)- The 
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FIG. 3. (Color online) Local density of states of Eq. (fldl near (a) the 
vortex and (b) the antivortex of the {dx 2 _y 2 + idxy)-wsL\Q pairing 
superconductor. The distance from the (anti)vortex core is denoted 
by r. The intrinsic line width is taken to be T = 0.002. 


(a) 


Vortex 
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FIG. 4. (Color online) Local density of states of Eq. d near (a) 
the vortex and (b) the antivortex of the {dxz idyz)-wsL\Q pairing 
superconductor. The intrinsic broadening is taken to be T = 0.002. 


energy spectrum of the vortex-bound sates is found to be (see 
Appendix IaU^ 





( 12 ) 


where I = with k G N. Hence, in agreement with the 

topological argument of Refs. [H d, [ssl andlH, we find that 
the {dx 2 _y 2 + idxy)-v^SiyQ superconductor does not have any 
zero-energy vortex-bound states [Fig. [Ha)]. The first subgap 
state has non-zero energy e i and the other low-lying states are 
evenly spaced with spacing ei. 

We observe that the bound-state energy spectrum of the vor¬ 
tex is the same as the one of the anti-vortex. But there is a 
striking difference in the wavefunctions between the vortex- 
and anti vortex-bound states. This is demonstrated in Fig. [2l 
which plots the wavefunction amplitude P(r) = \u{r)\‘^ -k 
|t;(r)p of the first two lowest energy vortex- and antivortex- 
bound states. Indeed, from the above discussion and using the 
fact that (g^r)]*, we find that the wavefunc¬ 

tion amplitude of the k-th lowest energy (anti)vortex-bound 
state is given by 


P(r) = (13) 

f Re[/(r)Pj^^\(gr)]2-l-Re[^(r)P^i\(gr)]2 vortex 

antivortex. 

We observe that Re[Pa^^ (r)] exhibits a node at r = 0 for all 
a except for a = 0, in which case Re[PQ^^(0)] = 1. Hence, 
it follows that for the vortex-bound states the lowest-energy 
wavefunction (k = 1) is peaked at finite r, whereas for the 
antivortex it is peaked at the origin r = 0, see Figs. [2ta) 
and[2tb). This finding is corroborated by our numerical simu¬ 
lations, which we present in the following subsection. 


B. Numerical results 

To compute the energy spectrum, the wavefunction am¬ 
plitudes, and the local density of states of the vortex-bound 
states, we discretize Hamiltonian ([T]) in the presence of the 
vortex/antivortex pair @ on the tetragonal lattice with 200 x 
200 sites in the xy plane and 100 points along the kz direc¬ 
tion. The eigenenergies and eigenfunctions (i*) of 
Eq. O are obtained by exact diagonalization of the disrec- 
tized Hamiltonian. Since the system has translation invariance 
along z we can diagonalize it for each kz separately. In the 
following we fix the parameters to t = 1, p = 0.3, /x = 1.9, 
and p = 5. For the {dx 2 -y 2 + X(ia^^)-wave pairing state we 
chose Ao = 0.05, while for the {dxz + idyz)-y/ 2 iYe state we 
set Ao = 0.2. With this parameter choice, the Fermi surface 
has a cigar-like shape elongated along the kz direction, which 
corresponds to the shape of the Fermi surface at the H point 
of SrPtAs2i“i2i. We have checked that different parameter val¬ 
ues do not qualitatively change the energy spectrum and the 
local density of states of the (anti)vortex-bound states. 


1. Energy spectrum and wavefunction amplitude 

In Fig. [T] we present the energy spectrum Sy^k^ a function 
of kz of the {dx 2 -y 2 -h X(ia^^)-wave and {dxz + idyz)-y/^yQ 
pairing states. The energies of the vortex-bound states are 
indicated in red, showing that vortices in the {dxz + '^dyz)- 
wave state exhibit a zero-energy fiat band of bound states, 
whereas vortices in the {dx 2 -y 2 + X(ia^^)-wave state support 
bound states only at finite energy. 

Fig. [2] displays the wavefunction amplitude | (r) | of 

the lowest and second lowest energy bound state at the vor¬ 
tex and antivortex of the {dx 2 -y 2 + idxy)-y^ 2 iYe superconduc- 
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FIG. 5. (Color online) Broadened local density of states (LDOS) near 
the cores of (a) the vortex and (b) the antivortex of the {d^ 2 _y 2 + 
idxy )-^pairing superconductor. To simulate the effects of finite 
temperature and finite experimental resolution the LDOS was convo¬ 
luted with a Gaussian with FWHM r = 0.045 corresponding to the 
energy resolution of the experiment. 


FIG. 6. (Color online) Broadened local density of states (LDOS) near 
the core of (a) the vortex and (b) the antivortex of the {dxz + idyz)- 
wave pairing state. The LDOS was convoluted with a Gaussian with 
FWHM T = 0.061 corresponding to the energy resolution of the 
experiment. 


tor for kz = 0 . Our numerical results (dashed blue curves) 
are in excellent agreement with the analytical solutions of the 
(anti)vortex-bound states (solid red curves). As discussed in 
the previous subsection, the lowest energy antivortex state ex¬ 
hibits a peak at r = 0, whereas the vortex state peaks at a 
non-zero r. For the second lowest energy bound state the be¬ 
havior is opposite: The vortex state is peaked a the origin, 
while the maximum of the antivortex state is at finite r. We 
note that these trends are independent on the kz value, since 
the overall shape of the wavefunctions is given by the order 
k of the Hankel functions which only depends on the 
vorticity n and the quantum number I (see Sec. lIII Al) . 

We remark that there is a similar asymmetry between the 
vortex- and antivortex-bound states of the {dxz + idyz)-^^^^ 
pairing superconductor. That is, the zero-energy anti vortex- 
bound state has a maximum at r = 0, while the zero-energy 
vortex-bound state is peaked at finite r. Again, this is a conse¬ 
quence of the difference in the order k of the Hankel functions 
describing the bound-states of the (anti)vortex (cf. dis¬ 
cussion in Refs. [13-113) • 


2. Local density of states 

The vortex-bound states of type-II superconductors can be 
probed by scanning tunneling spectroscopy of the surface den¬ 
sity of states^"—. To facilitate direct comparison with experi¬ 
mental measurements, we calculate the local density of states 
(LDOS) around the vortex and antivortex cores. The local 
density of states as a function of distance r from the vortex 


(or antivortex) center is given by 





E iT — 


(14) 


where v labels the eigenstates and eigenvalues e^^kz ’ ^ 
denotes the energy, and F represents an intrinsic broadening 
due to disorder. 

Figures [3] and [4] show the LDOS near the core of the vor¬ 
tex/antivortex of the {dx 2 -y 2 -h idxy)-w 2 iYe and the {dxz + 
idyz)-wsiyc states, respectively. The bound states appear as 
sharp peaks with energy spacing ei, given by Eq. (O. Com¬ 
paring Fig. |2| with Figs. [3] and m we find that the non-zero kz 
dispersion of the finite-energy bound states leads to a small 
broadening in energy of the LDOS peaks. The Majorana flat- 
band states of the {dxz + idyz)-w 2 iYe paring superconductor, 
on the other hand, have no kz dispersion and therefore give 
rise to very sharp zero-energy peaks near the center of the 
vortex and antivortex cores, see Fig.[4l Importantly, the asym¬ 
metry between the vortex and antivortex-bound states is di¬ 
rectly visible in the local-density of states: The lowest-energy 
antivortex-bound states are peaked at r = 0, while the lowest- 
energy bound states of the vortex have nodes at r = 0. The 
reason for this distinction was discussed in Sec. lIII Al 

In a scanning tunneling scpectrocopy experiment the LDOS 
is smeared by temperature broadening^—. To simulate this 
we convolute the LDOS with a Gaussian with full width at 
half-maximum r corresponding to the experimental energy 
resolution. We choose r to be of the order of two times the 
level spacing of the bound states ei, Eq. (O. Figures [5] and [6| 
show the broadened LDOS for the two pairing symmetries. 
We observe that the LDOS peak of the vortex is much broader 
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than that of the anti vortex. Moreover, the peak of the vortex 
is about half the height of that of the antivortex and it exhibits 
two ridges which disperse away to larger r. This is because 
the height of the LDOS peak is determined by the broadening 
of the lowest-energy state, while for the vortex it is due to the 
broadening of several low-energy states. In conclusion, we 
find that the asymmetry between the vortex and the antivortex 
can be detected in the LDOS even at temperatures T larger 
than the level spacing ei. 


IV. DISCUSSION AND FINAL REMARKS 


In this paper we have used large-scale exact diagonalization 
and analytical methods to study the structure of vortex-bound 
states in chiral d-wave superconductors. We have shown that 
vortices in the chiral {dxz i idyz)-^diNQ state bind dispersion¬ 
less zero-energy states, which form a doubly degenerate Ma- 
jorana fiat band (Fig.[T]). The stability and robustness of these 
zero-energy vortex-bound states is guaranteed by a Chern- 
Simons topological invariant. For the dz idxy)-w2iYe 

superconductor we found that vortex-bound states exist only 
at finite energy. We have computed the LDOS near the core 
of the vortex and antivortex of these chiral d-wave supercon¬ 
ductors. Importantly, we found a pronounced asymmetry in 
the LDOS between the vortex and the antivortex: The lowest- 
energy peak in the LDOS of the antivortex has its maximum at 
r = 0 , while the lowest-energy peak of the vortex is centered 
at r 7 ^ 0 (Figs. [3] and HI). Moreover, we have shown that the 
Majorana vortex-bound states of the (dxz i idyz)-wsiye super¬ 
conductor give rise to a particularly sharp peak in the LDOS, 
since these zero-energy bound states do not exhibit any dis¬ 
persion in energy (Fig. SI). 

The asymmetry in the LDOS between the vortex and the 
antivortex can in principle be used as a clear experimental fin¬ 
gerprint of the chiral order parameter symmetry. In practice, 
however, this might naively only be possible at temperatures 
T smaller than the level spacing ei between the bound states, 
since the LDOS is smeared by temperature broadening. The 
energy spacing ei [cf. Eq. ([T^ l is of the order of Aq/Ef, 
where Aq is the superconducting gap amplitude and Epis the 
Fermi energy. For a typical unconventional superconductor 
this corresponds to a temperature of about ^ 100/iK, which 
is below the reachable temperature regime of current state- 
of-the-art STM machines^. However, a clear asymmetry be¬ 
tween the LDOS of the vortex and the antivortex remains even 
at temperatures T of the order of ei < T < Aq, see Figs. [5] 
and[ 6 l That is, even though the individual LDOS peaks of the 
bound states cannot be resolved at a temperature T > ei, the 
broadened peak around the vortex is much wider and about 
half as high as the one of the antivortex. Hence, we believe 
that the predicted asymmetry is experimentally accessible for 
realistic materials, such as URu 2 Si 2 ^“— and SrPtAs 2 i~i 2 Z, and 
hope that our findings will stimulate future STM experiments 
on these interesting unconventional superconductors. 
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Appendix A: Derivation of the vortex-bound states 


In this appendix, we derive analytical formulas for the so¬ 
lutions to the BdG equations ([TOb in the limit 7 ^ 1 . For 
brevity, we restrict our discussion to the ansatz 

u{r) = v{r)= gi{r)H^^\{qr) (Al) 


for the wavefunctions. The solutions for the ansatz in terms 
of the Hankel functions of the second kind can be de¬ 
rived in an analogous manner [cf. discussion above Eq. (IA9M . 
Assuming r ^ 1 /can approximate the Hankel function 


H^aXdr) 


exp[i(g'r — a7r/2)] 


(A2) 


We observe that the asymptotic form of is proportional 

to times a phase factor, i.e., = (—since 

M_ — M+ = —21 and I = In addition, we find that the 

derivatives of the asymptotic Hankel function (IA2I) are given 
by 


dH 


( 1 ) 

M+ 


dr 


d^H 


( 1 ) 

M+ 


dr‘^ 



(A3) 

(A4) 


Inserting ansatz (lAlb into Eqs. (fTOl) and using the above ap¬ 
proximations yields the following differential equations for fi 
and gi in the limit qr 1 


'2^m+ - 7 


h + 


-Lm_ + 7 


gi + 


7 

^(_l)k+l 

7 


D-h 


with the differential operators 


c/i, (A5a) 
egi, (A5b) 


L± = E 2iqdr 


4^2 


[5 + 4(M±)2+4gV2] (A5c) 


and 

(y-ig) (A5d) 

-\-h'^q ^ ^ ^ _ iqj.'j _ fi'fj 

The set of equations (IA5I) can be analyzed in a perturbative 
approach. In the small q limit and focusing on solutions that 
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are decaying as r ^ oc, we find that at the first order in q the 
equations are solved by the exponential functions 


/i(r) = exp ^[ao(r) + i/3o(r)] + ^[ai(r) + ,(A6a) 

gi{r) = exp ^[ao(r) - i/3o(?’)] + ^Wi{r) - /3i(r)]^ ,(A6b) 

with 

pT 1 

ao = — \/2h^{r')dr and /3o = —. (A6c) 

Jo 2 


The functions ai and in Eqs. (IA6b describe corrections at 
the next order in q and can be determined by approximating 
/i and gi by 


/i (r) « (^1 + J (ai (r) + (r)e-2«°) ^ , (A7) 

51 (r) « (-l)‘‘+ie“° (^1 + J (ai(r) - /3i(r)e-2«o)^ , 


equations for ai and /3i, we obtain 


= — f 3h^{r') — V2h{r')h\r')dr', 

Jo 

^i(r) = - ^ l^e - j (A8) 


We observe that the solutions /i and Eq. (IA6I) . with o^i 
and Pi given by Eq. (IA8I) . are well behaved for large r since 
the radial vortex profile h{r) approaches 1 at large distances. 

The coefficients /2 and g 2 for the Hankel functions of the 
second kind Ha\x) in Eq. iTT\i can be derived in a sim¬ 
ilar manner, repeating the same steps as above. We find 
/ 2 (^) = /i (^) and ^ 2 (^) = Einally, we are ready 

to construct the full solution to the differential equations ([TOb . 
which is given in terms of a superposition of and 
The full solution needs to be regular at the origin r = 0, which 
leads to the condition that lm[/i(0)] = lm[^i(0)] = 0. That 
is, /i(0) and gi{0) need to be the same for the two Hankel 
functions, such that the imaginary singular part of the Hankel 
function is eliminated at the origin. Erom Eq. (IA7I) we find 
that this requirement is equivalent to ai(0) = Pi{0) = 0. 
The condition for ai is automatically satisfied; the one for pi, 
however, yields 


I 


oo 



g2ao(r)^^ = 0, 


(A9) 


and substituting this ansatz into Eq. (IA51) . Equating terms which determines the energy spectrum of the vortex-bound 
which are q-independent and solving the resulting differential sates ei, which is given in Eq. 
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